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The unknown neutrino parameters may leave detectable signatures in the supernova (SN) neutrino 
flux. However, even the contribution from the MSW flavor transition alone could cause ambiguity 
in the interpretation to the neutrino signals because of the uncertain local density profile of the SN 
matter and the model-dependent SN neutrino spectral parameters. A specific parametrization to 
the unknown local density profile is proposed in this work, and the contribution from the standard 
MSW effect is investigated through a multi-detector analysis of the SN neutrinos. In establishing 
the model-independent scheme, results based on the existing spectral models are included. The 
limitation of the analysis is also discussed. 

PACS numbers: 14.60.Pq, 13.15.+g, 97.60.Bw 
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As a distinctive type of the neutrino source, the core-collapse supernova (SN) provides a rich 
physical content that is lacking in the terrestrial environment. With its unique production and detection 
processes, the neutrino burst from a SN has long been considered as one of the promising tools for 
probing the unknown neutrino intrinsic parameters [H, S S 0; H H 0| i m particular, the neutrino mass 
hierarchy and the tiny mixing angle #13. However, difficulties also arise from the complexity caused 
i-C ' by the unavoidable astrophysical uncertainties, which could lead to ambiguous interpretations of the 

observed events. 
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INTRODUCTION 



The existing paradigm for the neutrino MSW flavor conversion in a SN has been established with 
the assumption of small neutrino self interaction. However, a possible new paradigm has been shaped in 
recent years. The argument is that, in addition to the standard MSW effect, the neutrinos may encounter 
certain non-MSW effects in the SN environment when the neutrino number density is extremely high. 
These not so well-known effects may be independent of the MSW effect and may introduce additional 
factors that alter the efficiency of the neutrino flavor transition. The current consensus is that both 
types of effects could contribute to the flavor transition and should all be included in the more complete 
t^J" ■ analysis of the SN neutrino signals. 

CV 

Even the standard treatment of the MSW flavor conversion alone is not immune from ambiguity. 
(N ■ The uncertain local density profile of the SN matter and the undetermined neutrino spectral parameters, 

such as the average energy and the luminosity of each neutrino flavor, are the crucial factors involved 
in the MSW effect of the SN neutrinos. However, an analysis based on the simple global density model: 
p ~ r~ 3 , which is widely adopted in the literature, could lose the generality if the variation of the local 
density shape near the resonance is not taken into account. With these uncertain factors, it is worth 
while to investigate how the contribution from the standard MSW effect should be modified in analyzing 
the SN neutrino signals, unless the overall contribution from the non-MSW effects is much greater than 
that from the MSW effect. 

The resonant neutrino flavor transition in a SN and in Earth would in principle give rise to observ- 
able signatures that reflect the properties of neutrino parameters. The unknown local density profile 
of SN matter near the resonance takes part in the adiabaticity parameter of the level crossing, and 
plays a role in the determination of transition probabilities. In fact, the adiabaticity of the level crossing 
could vary abruptly with the local density profile in certain neutrino parameter space. In addition, the 
knowledge to the primary spectrum for each neutrino flavor is essential in assessing the efficiency of 
neutrino flavor conversion. Various mixing scenarios, which arise from the uncertain SN physics, the 
possible neutrino mass hierarchies, and the uncertain magnitude of #13 lead to distinct neutrino survival 
probabilities. An analysis of the observed neutrino signals should, in principle, be able to single out the 
working scenario. 

The promising features of the multi-detector experiments have been well recognized IqL Hoi. [Til. [l2| ■ 
With the uncertain local density profile and the spectral parameters, this work investigates the potential 
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signatures that may be related to the contribution from the MSW effect in a multi-detector analysis of 
the SN neutrinos. As suggested by the previous work [HI, [HI, the consequences due to the uncertain 
local density profile near the resonance can be accounted for by the adoption of independent and variable 
power-law density functions p(r) = Ckr nk . However, the calculations in Ref. 13 and 14 are performed 
under the assumption of nearly constant Ck , which is similar to that adopted in the usual analysis based 
on p{r) = cr~ 3 . In the present work, the key improvement in dealing with the variable density profile is 
the establishment of variable Ck , which reproduces reasonable densities that agree with the numerical 
simulations at a wide range of the radial location. With the more general parametrization of the density 
function, this work is further devoted to analyzing the expected neutrino events at two water Cherenkov 
detectors. Focus is aimed at the modified contribution from the MSW effect. Certain physical observables 
derived from the expected event rates at the two detectors are proposed as the discriminators of various 
transition scenarios for the MSW resonance. In searching for the model-independent properties of the 
observables, three existing SN neutrino spectral models proposed by the Garching group fl5l [l6j and 
the Lawrence Livermore group [ljj are adopted in the calculation. 



II. UNCERTAIN DENSITY PROFILE OF SN MATTER 



There is no definite way of modeling the practically unknown local variation of density profile. 
Considering the narrow thickness of the resonant layer, as compared to that of the whole scope of the 
SN matter distribution, the density profile with a form of variable power-law in each resonant layer would 
be a reasonable simplification that accounts for the possible dynamical consequences. The essence of 
this approximation is that the power characterizing the density profile in each resonant layer is allowed 
to vary independently: p(r) ~ Qr" ! for the lower resonance layer, and p(r) ~ c/ l r" h for the higher 
resonance layer, where the factor Ck denotes the magnitude of the profile r" fc , with k = I or h. 

In this present work, the density profile in the resonance layer is parameterized as 

,M = [^r, (i) 

where c is a mass scale, and i?o is a distance scale, which is conveniently taken as the solar radius, 
Rq ~ 6.96 x 10 10 cm. With the variable rik, an important criterion for a reliable parametrization of the 
density profile is that the predicted density at a specific radial location agrees with that given by the 
numerical simulation, at least in order of magnitude. This suggests that in Eq.(l) the variation of Uk is 
accompanied by the variation of c at a different location. One may thus write the density profile in a 
more convenient form as 

P{r) = [jk][y][^- } nk , (2) 

where cq is a constant mass scale, and cq/Rq is a reference density scale. In the typical model with a 
fixed power, p(r) — cr~ 3 , the magnitude of c varies weakly fl8j: 10 31 g < c < 15 x 10 31 <7 in the density 
range 10 _5 g/cm 3 < p < 10 12 g/cm 3 . For variable rik near a given location, however, the mass scale c 
may vary by several orders of magnitude so as to reproduce reasonable density at this specific radial 
location. Note that with this variable parametrization of the density profile, the choice of cq does not 
alter the results since the variation of the mass scale is represented by c. In this present analysis we 
shall adopt the value cq — 7.0 x 10 31 <?, which is simply the mean value of c in the typical density model. 
The variables c and r may be rewritten as c/cq = X and r/Ro = Z, respectively, and the density profile 
becomes 

p{r) = [^]XZ n K (3) 

By fitting the predicted densities with the numerical results 0,113, HH at different radial locations, this 
parametrization leads to an approximate relation that regulates the variation of c, r, and nj,: 

X~ ^Z-<< nk+3 \ (4) 

As summarized in the following, the variation of nk can lead to non-trivial effects on the crossing 
probabilities at the resonance. One may first write the electron number density N e (r) for a typical core 
collapse SN as 

N e (r) = p(r)[Y e /m n ], (5) 
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where the electron number per baryon Y e = 1/2 is adopted, and m n is the baryon mass. The adiabaticity 
parameter for the resonant transition, defined as[22j 



26i 3 



lk -2Ecos29 l3 \jr^W ^ 



becomes 



1 Sm 2 sin 2 20a, r/r , , , cos20« i (1+^-) , v 

iy = m 2 — to 2 , 0y is the mixing angle between the eigenstates Vi and i/j, Gf is the Fermi 
constant, E is the neutrino energy, and Z$ — ro/Ro indicates the location of resonance. Note that 
\(l/N e )(dN e /dr)\o is evaluated at the resonance, where 

<5m?- cos 26 a 

(N e ) = — ^ (8) 

V ' 2V2G F E V ; 

The adiabaticity parameter jf. appears in the level crossing probability P k (with k = I or h) as: 

exp[-f lk F k ] - exp[-f 7*3^77] 



l-exp[-f7 fe ; 



(9) 



where F k is the correction factor to a non-linear profile. The origin of Eq. (9) and the correction function 
F k can be found in, e.g., Ref. 22. We shall adopt 5m 2 21 = 7.0 x 1(T 5 eV 2 , |<5to 2 3 | = 3.0 x 10~ 3 eV 2 , and 
sin 2 26*12 = 0.81 in the calculation. With a given radial location for the resonance, it can be verified that 
the adiabaticity parameter j k , arid thus the neutrino survival probability, vary weakly with the energy 
E, as compared to the influences from the variations of n k and #13 . This result suggests that one may 
simply adopt the average neutrino energies, e.g., E Ve — 12 MeV and E^ c = 15 MeV, in the following 
calculation for the neutrino transition probability in SN. The energy dependence would only appear in 
calculating the Earth effect. 



III. SURVIVAL PROBABILITIES AND SCENARIOS FOR NEUTRINO PARAMETERS 

When the neutrinos arrive at Earth, the survival probabilities for v e and v e are given respectively 

by 

Pnor = UZtftPh + Ul 2 {l - Pl)P h + C/ 2 3 (l - P h ), (10) 
Pnor = Ut l {l-P l ) + U 2 e2 P h (11) 

for the normal hierarchy, and 

P l nv = U 2 e2 {l-P l ) + U 2 el P h (12) 



Pinv = U^PtPh + U 2 el (l - Pl)P h + U 2 3 (l - Phi (13) 

for the inverted hierarchy, where Ph(Ph) and Pi (Pi) represent the higher and the lower level crossing 
probabilities for u e (v e ), respectively, and U e t is the element of the mixing matrix. 

In calculating the above level crossing probabilities, it should be emphasized that both the higher 
and the lower level crossings occur in the v sector for the normal hierarchy, while the higher crossing 
occurs in the v sector and the lower crossing occurs in the v sector if the mass hierarchy is inverted. 
One thus needs to calculate Ph and Pi in Eq.(10) for the normal hierarchy, and Ph in Eq. (13) and Pi 
in Eq.(12) for the inverted hierarchy at the individual resonance. Note that since the antineutrinos do 
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FIG. 1: The rapid variation of Ph as functions of sin 2 20\-j, for rih = —8 and rih = —2. 




TABLE I: The survival probabilities for u e and v e , the possible neutrino mass hierarchy, and the function 
g{iih,6\3) as predicted by the specific scenarios. The uncertainty in the power of the local density profile can 
alter the predicted bound of #13 through the function g(rih, #13). 
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not encounter the resonance at the location where the neutrinos cross the lower resonance, the crossing 
probability Pi for both hierarchies in Eq.(ll) and Eq.(13) is vanishing: Pi = 0. In addition, it can be 
verified that the adiabaticity parameter 7/ at the lower resonance is very large, 7; > 10 2 . This suggests 
that the lower resonance is adiabatic for both the mass hierarchies: ji ^> 1, which leads to Pi ~ 0. 

On the other hand, to show that the variation of rih and #13 leads to non-trivial structures of Ph 
in the rih — 6*13 parameter space, we first estimate the typical density near the location of the higher 
resonance. Using Eqs.(5) and (8), the density at resonance is given by 



l§! COS 26»i3 



[2V2G F E][Y e /r 



(14) 



which leads to p res ~ 3x 10 3 g/cm 3 with small #13. This density corresponds to an approximate location 
of resonance tq/Rq = Zq ~ 0.03 according to the numerical results [19|, [20j, |2l| . It follows from Eqs.(7) 
and (9) that the crossing probability Ph is in general non- adiabatic (P h ~ 1) for sin 2 26»i 3 < 10~ 5 . This 
leads to P nor = F m „ = |t/ e 2| 2 ~ 0.3, as can be checked using Eqs. (10) and (12). In addition, it can 
be verified that Pi nv remains constant through out the parameter space with Pi nv ~ 0.3. However, the 
probability P nor varies with both rih and #13 near sin 2 2#i3 ~ 10~ 5 — 10~ 4 since Ph drops rapidly from 
Ph ~ 1 (non-adiabatic) to Ph ~ (adiabatic) within a narrow region of the parameter space. It follows 
that Pnor drops from P nor ~ 0.3 to P nor — \U e 3\ 2 <C 1. This property distinguishes P n0 r from Pi nv . As 
an illustration, we show Ph as a function of sin 2 20i3 for nh = — 2 and rih — — 8 in Fig. 1. 

It is then reasonable to represent the narrow borderline between the adiabatic crossing (Ph ~ 0) 
and the non-adiabatic crossing (Ph ~ 1) by the condition Ph = 1/2, which implies (with small #13) 

exp[-~ lh F h ] ~ I. (15) 

Equations (7) and (15) lead to jW2^ h = 9( n h,6i3) ~ lj with 

7T 1 |jm|i I sin 2 26<i3 
9{nh,V 13 ) = - — 1 — — 

4(ln2) \rih\ & cosiuis 

xlQSmll/E) ;° S l 9 :; — ]^R Z^\ (16) 

where Zq — 0.03, Ro = 6.96 x 10 10 cm, and Fh ~ 1 (for small 6*13). Thus, how the uncertainty in rih 
alters the predicted bound for 6*13 is regulated by the function g(nh, O13), and depends on whether the 
parameters rih and 813 result in 3(71^,6*13) < 1 or g(nh,0\?,) > 1. The similar properties for P nor and 
Pinv also can be derived from Eqs. (7), (9), (11), and (13). As a brief summary, one notes that with 
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FIG. 2: The adiabaticity parameters is modified for a finite resonance layer, within which the density and the 
mass scale c are both variable. The contours for g(nh, #13) = 2 ™ 2 fh = 1 are shown in the nu — #13 space with 
Z = Z , Z = Z (l + 10%), and Z = Z (l - 10%), where Z = 0.03. Note that since the estimated deviation 
(< 2.7%) is barely observable in the parameter space, an exaggerated uncertainty of 10% is shown here for the 
purpose of illustration and comparison. We may conclude that this small deviation does not alter the general 
picture of the analysis. 




the formulation and the chosen input parameters, (i) Pi — 0, and Pi — 0, and (ii) Ph = Ph — 1 if 
g(n h , 013) < 1, while P h = P h = if g(n h , 9 13 ) > 1. 

In terms of the survival probabilities P and P for v e and v ei respectively, the uncertain local 
density profile, the possible mass hierarchies, and the undetermined mixing angle #13 together lead to 
three possible combinations of (P,P), as indicated by (a), (b), and (c) in Table HI In general, separating 
the scenarios is unlikely if sin 2 26 13 < 1CT 5 since P n0 r and P m „, as well as P n0 r and Pi nv , cannot 
be distinguished. In addition, there is a non-trivial dependence of P and P on both rih and #13 near 
sin 2 2013 ~ 10~ 4 . This implies that establishing a better lower bound for #13 relies on whether the 
uncertainty of nu in this region of parameter space can be reduced. Based on Table HI we discuss the 
following cases for the purpose of illustration, (i) If the mass hierarchy is identified as normal by the 
observation of extremely small P, the predicted lower bound for #13 could still be uncertain by more 
than one order of magnitude because of the uncertain nh- In this case, the condition g(rih,0i3) > 1 
establishes the predicted lower bound for #13 with a specific value of rih- (ii) On the other hand, a 
small P implies inverted hierarchy, and the predicted bound for #13 is also regulated by the condition 
g(nh, 613) > 1- (iii) If both P and P are observed to be large: P ~ 0.3 and P ~ 0.6, then it implies 
g{nh,6\3) < 1 with no information about the mass hierarchy. 

Note that with the estimated location of resonance Zq = 0.03, there is a potential source of un- 
certainty in calculating the adiabaticity parameters. The key point is, with the variations of nh, c, and 
the density in a finite resonance layer, what modification to 7^ is required? We first show in Fig. 2 
the contour (solid line) of g(rih, #13) = 1 in the n/, — #13 space for Zq = 0.03. In order to estimate the 
uncertainty, one needs to calculate the width of the resonance layer 

2tan26»i 3 

I JV e dr 1° 

which can be reduced to So ~ 2(0.03i?o) tan 2^13/ |nh|. One notes from the solid contour of Fig. 2 
that, as nh varies in the region —10 < n% < —1.5, the higher resonance occurs roughly in the region 
5 x 10~ 4 < sin 2 26<i3 < 4 x 10~ 3 . This leads to an estimated bound for the width, 5 /(Z R ) < 2.7%, 
for the given nh and $13 in the parameter space. The deviation of the contour g(rih,@i3) = 1 caused 
by this small width is barely observable in the nh — 6*13 space of Fig. 2. To illustrate the smallncss 
of this deviation and as a comparison, we also show in Fig. 2 the contours of g(rih, #13) = 1 due to 
an exaggerated width of the resonance layer, regardless of the values for any given rih and #13. It is 
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TABLE II: The average energies and relative luminosity for neutrinos of different flavors suggested by the Gl, 
G2, and the LL models. 
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FIG. 3: The neutrino incident angles at detectors A and B, as measured individually from the local zenith, are 
given by ip a and ipp, respectively. The angle formed by A, O (the Earth center), and B is denoted as A. 




seen that even with the conservative estimation based on the enlarged width Z = Zq(1 ± 10%), the 
resultant deviation of the contour from that of Zq = 0.03 is still very small. We may conclude that the 
small uncertainty due to the finite width of the layer does not alter the general picture of the analysis. 
Note that all the three contours meet at Uh = —3. This can be realized from Eq. (16), in which the 
contribution from any deviation of Zq vanishes when rih = —3. 

The survival probabilities will be modified by the regeneration effect as the neutrinos propagate 
through the Earth matter. The three possible combinations (a), (b), and (c) for (P, P) are now modified 
according to 



P( Q ) = sin 2 13 P 2e , 

P (b) = P 2e , 

P (c) = P 2e , (18) 

and 

P< Q ) = 1-P 2e , 

PW = sin 2 Ml-P 2e ), 

P^ = 1-P 2e , (19) 

where P 2e (P 2e ) is the probability that a z^ 2 (P2) arriving at the Earth surface is eventually detected as 
a v e (pe) at the detector [H, BH- 



IV. SN NEUTRINO SPECTRAL PARAMETERS AND MULTI-DETECTOR ANALYSIS 



In addition to the brief pulse of the neutronization v e burst, all three flavors of neutrinos and 
antineutrinos are emitted from the SN through pair production processes during a typical time scale 
~ 10s. The primary SN neutrino spectrum, which is not pure thermal, is usually modeled as a pinched 
Fermi-Dirac distribution with several spectral parameters: 



where Li is the luminosity of the neutrino flavor z^, 7j is the effective temperature of v% inside the 
respective neutrinosphere, E is the energy, feiyji) is the normalization factor, and r\i is the pinching 
parameter for F®. Notice that r} Ue ~ rjp e ~ r/^ x ~ 3.0 is assumed here, where u x = u^, v T , v^, and v r . 
The most significant difference between the Garching models (Gl, G2) and the Lawrence Livermore 
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FIG. 4: The ratio & a ' as a function of the neutrino energy, with ipa = 7r/4 (detector A is unshadowed) and 
t[>0 = 0.957T (detector B is shadowed by mantle and core). The primary neutrino spectra are adopted from Gl, 
G2, and LL models. 




model (LL) is that the Garching models exhibit a less hierarchical structure in the luminosity and in 
the average energy among different neutrino flavors, as shown in Table 1111 

We assume a core collapse SN at a distance 10 kpc away. In a typical water Cherenkov detector, 
the isotropical inverse /3-decay, v e + p — > n + e + , is the most dominant process. In this work, other 
isotropical CC processes, v e + O — > F + e~ and v e + O — > N + e + , are also included in the calculation. 
As for the directional events, we include the scattering events induced by Vi and v e : Vi{v e ) + e~ , where 
i = e,/i,r. The cross sections for P e +p and Vi(9 e ) are adopted from Ref . [2a] . and that for v e (T>e) + O 
are adopted from Ref.[26j]. As shown in Fig. 3, the neutrino incident angle at a detector, denoted as tp a 
and for detector A and B, respectively, may be defined as the angle measured from the local zenith. 
For a given i/j a , the allowed range for ipp is limited. If the angle formed by A, the Earth center, and B 
is given by A, then 

|V>a - A| < Ipp < + A (21) 

for ip a + A < 7r, and 

|^a-A| < i>p <2ir-(^ a + \) (22) 

for 7r < tp a -\-X < 2tt. The density of Earth matter is of the order ps ~ lOg/cm 3 , and its variation is much 
smaller than that of the SN matter. One may adopt a simple step function [2 7j| as the approximation: 
Pe ~ 5.0 g/cm 3 for ir ffi < r < r® (mantle), and pe ~ 12.0 g/cm? for r < ^r® (core), where r® is 
the Earth radius. For the purpose of illustration, one may analyze the expected results at two future 
detectors: the HyperKamiokande [28j (detector A) and the MEMPHYS(29j (detector B). The locations 
of the two detectors lead to A ~ 7r/2. 

Since the neutrino flux penetrates significant amount of Earth matter only when the incident angle 
at a detector is greater than tt/2, no observable difference is expected from the two detectors if both 
%j} a and ipp are less than tt/2. In this case, both Pi e and p2 e reduce to sin 2 #12, and all the observed 
differences between A and B vanish. Given the variety of spectral models, the present analysis is aiming 
at seeking for model-independent properties that would help identify the working scenario. 



V. OBSERVABLES AND ANALYSES 



Simultaneous observation of the SN neutrinos at two terrestrial detectors would yield the most 
useful results when pronounced Earth regeneration effect is observed at only one of the detectors. 
Given the two terrestrial detectors A and B, a proper analysis of the expected neutrino event rates at 
the two detectors would in principle reveal the signatures representing distinct scenarios of the neutrino 
properties and mixing schemes. To pave the way for the analysis, we first examine the energy dependence 
of the relative flux difference at the detectors A and B: 
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for the v e flux, and 



pB _ pA 

D = VC F A V * ( 23 ) 



pB _ pA 

D = *'p_A ( 24 ) 



for the v e flux. With F®. and F®. denoting the primary flux of Vi and Vi, respectively, one may relate 
the observed and the primary fluxes as 



F Ve =Fl + {l-P){F° x -F° e ), (25) 



F De =Fl + {l-P){Fl-Fl). (26) 



The three scenarios then lead to several distinct relative flux differences: 



D(a) = sm 2 6 13 lF^-F°JlP 2e -sm 2 e 12 ] 

F° [sin 2 9 13 sin 2 12 ] + F° [f - sin 2 13 sin 2 12 ] ' 



[F° -F°][P 2e -sin 2 



F° sin 2 12 + F$ cos 2 12 ' 



D(c) = [F° -F°J[P 2e -sin 2 g 12 ] 
F° sin 2 12 +po cos 2 12 ' 



(27) 



D(6) = ^ ^ ^ (28) 



(29) 



£)(«) = 



-[^-^J[P2e-sin 2 g 12 ] 

F§ e cos 2 12 + Fl sin 2 12 : 



(30) 



-sin 2 13 [F° ' -F p °J[P 2e -sin 2 



F£ [sin 2 6»i3 cos 2 12 ] + F$ [f - sin 2 13 sin 2 12 ] ' 



(31) 



(c) = -[F° -F°J[P 2e -sin 2 12 ] 

F£cos 2 0i2+F£sin 2 12 ' 1 ' 

It is seen that the ratios also depend on the chosen model for the primary V{ and vi fluxes, as well as 
on P 2e and P 2e , which introduce the non-trivial energy dependence into D and D. As an illustration, we 
show the variation of as a function of energy in Fig. 4, with ip a = 7r/4 (detector A is unshadowed) 
and tpp — 0.957T (detector B is shadowed by mantle and core). The primary neutrino spectra are adopted 
from the Gl, G2, and the LL models. Note that with the finite resolution power of the detector, not all 
the details of the curve can be observed in practice. It is clear from Eq.(30) that the energy dependence 
of lM a ' originates from the factor P 2e — sin 2 0i 2 . The energy dependence of D and D, as shown in 
Eqs. (27-32), can be used to illustrate the properties of the observables that shall be proposed in this 
work. 



A. Relative difference of the event rates 



In terms of the event rates, one may first examine the relative difference of the observed event rates 
per unit target mass for detectors A and B: 

SS M. ,33, 
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FIG. 5: The possible range of the value for R = R(E) is indicated by the vertical range between two horizontal 
lines. Results for scenarios (a), (b), and (c) are shown here with the chosen bin size of 5 MeV. The incident 
angles ip a < n/2 (detector A unshadowed) and = 0.95-7T (detector B shadowed by mantle and core) are used 
to optimize the result. The estimated ranges of R shown in the figures are determined by the results of the three 
spectral models and the statistical errors. 
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The total event rates per unit target mass, Na and Nb, consist of all the isotropical and directional 
events mentioned in the previous section, with 

N A . B = (34) 

where Ma.b denotes the total target mass of a detector. The expected event rates 2Vj induced by the 
flavor Vi (or Vi) at the detector is given by 

at r 

N i = JzknE F »i*°iv]dE, (35) 
J v 

where N t is the target number at the detector, L is the distance to the supernova, Oi y is the cross section 
for Vi (or vi) in a particular reaction channel y. The detection efficiency is assumed to be one. In the 
following, one may denote the resultant values of R for the six scenarios as R J , where J = (a), (6), (c). 

As an illustration, we show in Fig. 5 the expected values of R for scenarios (a), (b), and (c). 
The reasonable estimation of the energy bin size is roughly a few MeV for the detection of typical 
SN neutrinos (30l|. We adopt an energy bin size of 5 MeV in the figure. Note that for scenario (6), 
the deviation of R from zero is ~ 1% or smaller, and the sign of R is undetermined with the error 
bar. The estimated width for R is given by the different results from the three spectral models and 
the statistical errors from each. As an example, one notes that the expected values of R^ under LL, 
Gl, and G2 models in 25 MeV < E <30 MeV are estimated to be -0.085 ± 0.057, -0.221 ± 0.112, 
and — 0.216 ± 0.109, respectively. The resultant R^ is then shown in Fig. 5 with the combined width 
-0.333 < i? (a) < -0.027 for 25 MeV < E <30 MeV. Note that there is no apparent energy dependence 
of the width. 

Even though some details of the energy-dependence in Fig. 5 are lost, certain useful qualitative 
properties are still available. It is seen that the Earth effect could suppress the event rates to an 
observable level (of order 10% or more) for E > 25 MeV if the working scenario is (a) or (c). The sign 
and the magnitude of R in E > 25 MeV can be used to separate the three scenarios into two groups: 
[(a), (c)] and (b). It should be pointed out that for E < 25 MeV, the uncertainties wipe out the small 
deviation from zero and the sign of R is undetermined for all scenarios. 

More hints may be available from analyzing the directional events alone, which are predominantly 
the t'e-induced elastic scattering events. However, the directional event rate is approximately 1^2 
orders of magnitude less than the isotropical event rate. Extracting the directional events effectively 
from the dominant background of the isotropical events may be quite challenging in practice. In the 
following analysis, one assumes that the directional event rates in the forward cone with a narrow solid 
angle can be identified. 
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FIG. 6: The values of S = S(E) estimated with the bin size of 5 MeV and ip a < n/2 (detector A unshadowed), 
ipp — 0.957T (detector B shadowed by mantle and core). The ranges of S shown in the figure are based on the 
combined results of the three spectral models and the statistical errors. 
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B. Relative difference of the directional events 



We may define another ratio based on the directional events per unit mass of target at the two 
detectors, 

n A 

where ub and tia denote the shadowed and the unshadowed directional events per unit target mass, 
respectively. The properties of S, as shown in Fig. 6, lead to new information that is unseen from R. One 
notes that two groups of scenarios, (a) and [(b), (c)], could be identified by the magnitude of S at the 
high energy end since the magnitude could differ by a factor of three or more. Thus, one may conclude 
that the degeneracy between and R^ could be removed by using the observable S, and that the 
three scenarios could be reasonably distinguished by analyzing both R and S at the high energy end. 

It should be pointed out that even though R and S represent energy-integrated quantities, they 
still exhibit general energy dependence with the chosen finite bin size, as can be seen from Figs. 5 and 
6. To qualitatively illustrate this general properties, one may approximate R, for simplicity, as the event 
ratios that consist only of the predominant D e +p events. For scenario (a), Eqs. (25-26) and (33-35) lead 
to 

„(„) „ [Fl -Fj]f(P 2e - sin 2 9 12 )a(E)dE 

j(Fl cos 2 6 12 + Fl sin 2 6 12 )a(E)dE' [ 1 

where cr(E) is the cross section. With the input parameters and a given model for the fluxes, 

(o) j(P 2e - sin 2 6 12 )a(E)dE 
R Ja(E)dE • (38) 

If one calculates R^ in a relative small step in energy, AE <~ a few MeV, the cross section a — cr(E), 
which is roughly a <~ E 2 , would be relatively smooth as compared to the function P 2e . As a consequence, 
one may expect ~ P 2e — sin 2 6>i2 . It is seen from Eq.(30) that the qualitative behavior of the factor, 
P 2e — sin 2 6*i2, which dictates the general trend of the energy dependence of in Fig. 4, propagates to 
energy dependence of S, and that of T, which is introduced in the next subsection, 
can be realized in a similar way. 



C. Double ratio with both directional and isotropical events 



The observables, R and S, tend to flip signs near the peak of the spectrum (~ 25 MeV), and become 
less useful for analyzing the abundant events induced by the neutrinos with energies near the peak. It 
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FIG. 7: The possible ranges of T = T(ipp) for scenarios (a), (6), and (c) are estimated from the combined results 
of the three spectral models and the statistical errors for 15 < E < 30 MeV. The incident angle for detector A 
is chosen as tp a = tt/2, and an assumed directional resolution of 7r/12 for detector B is adopted. The resultant 
range of T in each scenario is expected to show no -(/^-dependence for ipp < tt /2. 
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would be worth while to further investigate whether there is more to learn from the middle range of 
the spectrum. Theoretically, the Earth regeneration effect is signaled by the oscillatory modulations 
of the observed spectra. A quantitative analysis of the resultant event rates, however, may be quite 
subtle. The resolution power of the detector is limited and the modulation may not alter the total event 
rates significantly enough to render an observable difference between the shadowed and the unshadowed 
events. To optimize the results for an analysis, neutrino event rates from a specific range of energy cut 
may be needed. 

One may construct a double ratio with both the isotropical and the directional events: 

T _ [{Nb - n B ) - (Na - n A )]/(N A - n A ) ^ 
(n B -n A )/n A 

where N — n and n denote the isotropical and the directional events, respectively. Depending on the 
time of a day, a given incident angle ip a for detector A corresponds to a definite incident angle ipp for 
detector B. Notice that the possible values of ipp are regulated by Eqs. (21) and (22). If ip a ~ or 
ip a ~ 7r, the angle ipp is limited to a narrow range: ipp ~ A ~ 7r/2. On the other hand if ip a ~ 7r/2, the 
allowed range for ipp becomes maximized: < "0/3 < 7r, and the observables are expected to vary with 
a wide range of ijj/3- Each of the two detectors may be unshadowed (if ip a ,p < 7r /2), shadowed by the 
mantle (if tt/2 < ip a< /3 < 5ir/6), or shadowed by the mantle and the core (if 5ir/6 < ip a ,p < 7r)- In any 
case, more useful hints from a simultaneous observation would be available if one of the incident angles 
is greater than 7r/2 and the other is less that 7r/2. 

The observable T may be used as a supplementary tool to the other two observables, R and S. 
As an illustration, the ratio T is evaluated here for 15 < E < 30 MeV and analyzed as a function 
of ipp, with ip a = tt/2 and a conservatively assumed directional resolution of 7r/12 for detector B, as 
shown in Fig. 7. Notice that the pointing accuracy of the incident neutrinos at the SK detector can 
reach the level to within a few degrees at 10 kpc [3l|. The results for each scenario are expected to 
show no i/j^-dependence for ipp < tt/2. In addition, calculations show that scenarios (b) leads to small 
deviations from zero: —0.078 < T < 0.076, which can be distinguished from the other two scenarios. 
Furthermore, and can also be distinguished by the difference in magnitude. One would expect 
that a combined analysis of i?, S, and T should be able to provide certain useful constraint or prediction 
about the neutrino parameters. 



VI. SUMMARY AND CONCLUSIONS 



The SN neutrino burst has long been considered as a promising tool for probing the neutrino 
mass hierarchy and the mixing angle #13. The current consensus is that the contributions from both 
the MSW and the non-MSW effects should be included in analyzing the flavor transition of the SN 
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neutrino signals. However, the uncertain local density profile of the SN matter and the neutrino spectral 
parameters may lead to diverse interpretations to the neutrino signals from the contribution of the 
MSW effect. To investigate the consequences resulting from these uncertainties, this work suggests a 
more general parametrization for the uncertain local density profile, and outlines a scheme using ratios 
of experimental observables as the discriminators for the possible scenarios that are related to the MSW 
effect. As an illustration, we analyze the expected event rates at two terrestrial Cherenkov detectors. The 
consequences due to choices of the existing spectral models are included in the analysis. The possible 
scenarios arising from the uncertainties are then examined through the observables R, S, and T. It is 
seen that the uncertain local density profile of the SN matter can alter the predictions for #13 and the 
mass hierarchy from the contribution of the MSW effect alone. The predicted bound for 6*13 is regulated 
by the proposed function g{ri} ll 613). 

The results in this work are derived from the expected event rates at two megaton-scale Cherenkov 
detectors that would be available in the future. Analyses of the future SN neutrinos based on other 
types of detectors, such as the liquid Argon and kton-scale scintillation detectors that are sensitive to 
distinct channels of the neutrino interactions, should also provide valuable information toward a better 
understanding of the neutrino parameters and the SN physics. A detailed analysis based on all the 
three types of detector is performed in, e.g., Ref. 6, in which seven SN and neutrino parameters are 
taken into consideration. To extend the present work, it would be intriguing to investigate how the 
uncertain density profile would impact the outcome of the analysis based on the detectors other than 
the Cherenkov ones. 

In shaping a possible new paradigm for the neutrino flavor conversion in a SN, recent studies and 
simulations suggest that the shock wave propagation and certain types of non-MSW flavor conversion 
may occur in the complex environment of a SN at different space and time scales if the proper phy sical 



conditions are met. These effects, such as the neutrino self coupling in dense media[32L I33L 134 13 
37, 38 , [3^, H3, EH and the neutrino flavor de-polarization associated with the after shock turbulence 



43 



3G 
42 



44, 45]. may have different origins from that of the MSW effect and would impact the efficiency of 



neutrino flavor transition. 

It is suggested that the neutrino self interactions induce collective flavor transitions, which may 
occur before or in the same region as the resonant flavor transition. The size of modification to the 
original neutrino flux depends on the mass hierarchy, as well as on the mixing angle #13. As an example, 
one considers the simplest case when the self-interaction effects can be factored out from the ordinary 
MSW effects. For the inverted hierarchy, it is suggested that if #13 is nonzero, the collective pair- 
conversion of the type v e v e — > v x v x can be triggered before the ordinary MSW effect occurs[46|]. Thus, 
even in the simplest case, the altered primary spectra could impact the validity and the usefulness of the 
observables proposed in this work. A better understanding of the modified neutrino spectra resulting 
from the collective effects would help in establishing more convenient and useful observables since the 
knowledge of the neutrino spectra before and after the MSW effect is essential for the present analysis. 
Note that the treatment of the variable density profile for the MSW effect, as that proposed in this 
work, is unaffected by the non-MSW effects. 

The shock wave propagation, on the other hand, may alter the density profile of the SN matter 
and induce certain non-linear effects that modify the neutrino survival probabilities and the neutrino 
spectrum. In particular, there is a possibility that more than one MSW-resonance related to the same 
scale of mass-squared difference could arise due to the shock wave propagation. The effects of the 
shock wave on certain physical observables, with the emphasis on the time-dependent properties, are 
analyzed [H, [3(| S3|- Although the proposed scheme in the present work is limited to the discussion 
without considering specific scenario due to the shock effects, it takes into account the possible varying 
density profile, and is not limited to whether the profile is actually primary. In addition, the quantities 
i?, S, and T are proposed in this work to examine the time-integrated behaviors of the event rates. 
As an extension to this present analysis, it would be intriguing to further investigate weather the time 
structure of the shock waves effects and the possible multiple MSW-resonance due to the shock effect 
can alter the time-integrated properties of the proposed quantities to an observable level. Even with 
the unknown neutrino parameters, the uncertain model for the SN neutrino spectrum, and the possible 
statistical errors, the quantities R, S, and T might still be able to provide another potential means for 
probing the SN physics if they reveal time-independent signatures that are related to the shock wave 
propagation. On the other hand, if the shock signatures are unobservable with these time-integrated 
quantities, then these quantities might be able to provide insights to the neutrino parameters without 
being affected by the shock waves. 

It should be emphasized that the present work is devoted to effects that are related to the MSW 
resonant flavor transition in a SN. The focus is aimed at the possible modification to the results of 
MSW effect due to the uncertain density profile. This modification is expected to have certain impact 
to the complete analysis based on both the MSW effect and the non-MSW effects. The whole details 
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of the non-MSW effects, as well as the interplay between the MSW and the non-MSW consequences, 
are still not fully understood. The neutrino flavor conversion in a SN remains as a complicated problem 
that awaits for more complete answers. Whether the flavor transition arising from one origin would 
dominate over that from the others depends on the uncertain spectral parameters, the undetermined 
neutrino intrinsic properties, and the physical conditions that the neutrinos encounter at different stages 
of the SN environment. In addition, the reliability of the analysis relies on whether the MSW and non- 
MSW effects can be analyzed separately. However, these unrelated types of effects decouple only when 
specific physical conditions are met. 

In any case, analysis of the non-trivial MSW flavor conversion due to the possible drastic variation 
of the local density profile, as discussed in this work, should provide certain insight in probing the 
neutrino parameters and in shaping a more convincing paradigm for the neutrino flavor conversion in 
the SN. 
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